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Abstract 

Spider dragline silk shows well-known outstanding mechanical properties. However, its sigmoidal 
shape of the measured stress-strain curves {i.e. the yield) can not be described by classical polymer 
theories and recent hierarchical chain model. To solve the long lasting problem, we generalized the 
Maier-Saupe theory of nematics to construct an elastic model for the polypeptide chain network 
of the dragline silk. The comprehensive agreement between theory and experiments on the stress- 
strain curve strongly indicates the dragline silks to belong to liquid crystal elastomers. Especially, 
the remarkable yielding elasticity of the silk is understood for the first time as the force-induced 
isotropic-nematic phase transition of the chain network. Our theory also predicts a drop of the 
stress in supercontracted dragline silk, an early found effect of humidity on the mechanical property 
in many silks. 

PACS numbers: 81.05.Lg, 81.40.Jj, 82.35.Pq, 61.30.Dk 
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Spider dragline silks (SDSs), the main structural web silk regarded as the "spider's life- 
line" , exhibit fascinating mechanicalproperties, such as a tactful combination of high ten- 
sile strength and high extensibility [l|, thus showing a remarkably sigmoidal shape of the 
measured stress-strain curves Several experimental studies have been carried out to 
determine the supra-molecular structure organization of the SDS , [^l, Q] and tried to 

produce mimic silks with similar properties 8|. It is now widely accepted that SDSs are 
semicrystalline polymers with /3-sheet nanocrystals embedded in amorphous region, which 
is a polypeptide chain network [si, [lO, 11]; see Fig. 1(a). However, the deformation mech- 



anism, which is essential for understanding the SDS's extraodinary mechanical properties 



and mimicking the silk, is still in intense debate 



12 



HQ. 



On the theoretical side, to understand the exceptional mechanical properties of SDS is 
of longstanding interest, and many models have been proposed 0, 



insights attained 



flfl 
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12] . For example, the model by Termonia 9| treated the SDS 



111. Il2j| and some 



as a hydrogen-bonded amorphous region embedded with stiff crystals as cross-links. In the 
interfacial region, an extremely high modulus is required to get the dragline's overall behavior 



on deformation. While, in the model of Porter and Vallrath [Hi, |12|, parameters linking to 
chemical compositions and morphological order were used to interpret thermo-mechanical 
properties. But some parameters such as ordered/disordered fractions are difficult to be 



obtained from experiments. A recent model 



4j connecting deformations on macroscopic and 



molecular length scales still did not consider the change of the orientation of nanocomposites 

n 

during deformation. Especially, as pointed out by Vehoff at al. recently [13[, basic polymer 
theories such as the freely jointed chain, the freely rotating chain and the worm-like chain, as 
well as a hierarchical chain model of spider capture silk 15[ can not reproduce the sigmoidal 
shape or even the steep initial regime of the spider dragline silk [ Fig. 2(a)] 13]. In one 
word, a unified description for SDS as a model biomaterial still seems to be lacking. 

Quite a few works have pointed out that spider silk is liquid 

crystalline material and liquid crystal (LC) phase plays a vital role in both its spinning 
process and mechanical properties. In the spinning process, the liquid crystalline 'spin- 
ning dope' helps spider to control the folding and crystallization of the main protein con- 
stituents at benign condition (close to ambient temperatures and pressures using water as 



solvent) Q, 17]. The liquid crystalline phase also plays an important role in the solid 
silk's properties 181 ] . For instance, several works have found out that the orientation of 
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nanocomposites can affect SDS's mechanical properties significantly A recent ex- 

periment also suggested the existence of conformational transition and the liquid crystalline 
state of regenerated silk fibroin in water 8j. Therefore, to present an analytically tractable 
LC model of the SDS that can catch the main physical factors is a current challenge to the- 
orists. In this work, we generalized the Maier-Saupe theory [jjj] of nematic LC to construct 
an elastic model for the polypeptide chain network of the SDS. We show that on deformation 
the SDS undergoes significant changes with orientation of the chain network increased and 
the dimension of the silk along force direction elongated. The comprehensive agreement 
between theory and experiments on the stress-strain curve strongly indicates the SDSs to 
belong to LC elastomers, described as a new class of matter recently [2^. Especially, the 
remarkable yielding elasticity of the SDS is understood for the first time as the force-induced 
isotropic-nematic phase transition of the chain network and the self-consistently obtained 
yield point agrees with experimental data well. The present theory also predicts a drop 
of the stress in supercontracted SDS, an early found effect of humidity on the mechanical 



properties in many silks JjJ, 



13. 



2l|. 



We take the polypeptide chain network in the amorphous region of the SDS as a molecular 
LC field with each chain section corresponding to a mesogenic molecule; see Fig. 1. Because 
the SDS's high extensibility results primarily from the disordered region , [6[] , and many 
experiments showed that the deformation of the crystals is at least a factor of 10 smaller than 
that of the bulk and that the orientation of the /3-sheets is almost unchanged (usually 
very high) under stress lj], [3] , we can neglect the deformations and rotations of the /9-sheet 
crystals in current work. 

Following the LC continuum theory in the absence of forces, the potential of a mesogenic 



molecule takes the Maier-Saupe interaction form 

V(cos6) -- 



19] 



aS(^cos 2 8- -), 



(1) 



where 9 is the angle between the long axis of the molecule and the silk axis (the z-axis), 
which is also the direction of n [ Fig. 1 (b)], a is the strength of the mean field, and S 
is the orientation order parameter of the LC, defined as the average of second Legendre 
polynomial 23] 



S 



3 2n 1 

- cos 6 

2 2 



(2) 
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We notice that the Maier-Saupe potential has been used by Pincus and de Gennes in in 



vestigating LC phase transition in a polypeptide system 



24J. When a uniform force field f 



along z-axis is applied, the potential of a molecule is written as 

f/(cos0) = V-flcos6 : (3) 

where / denotes the length of the mesogenic molecule. 

From the definition of the order parameter S, we get a self-consistency equation 

S = /_!(§ cos 2 8 - ~) exp(J^ cos 2 8 + acos9)dcos8 
I § X _ X exp(J^ cos 2 8 + a cos 8)d cos 8 , 

with a = fl/kftT. The solution of the above equation may not be unique, in order to obtain 
physically sound solution we still need the requirement of minimization of the free energy 
given by 

Fms = -k B T In Z + ^aS 2 , (5) 

where Z is the partition function Z = f\ e~ u ( cosd ^ kBT dcos8, and the second term at the 

□ 

right-hand side corrects for the double counting arising from the mean field method [201 ] . 

We calculate the orientation function S numerically at temperature T* = T/T n ; and 
force / and show the results in Fig. 2(b). Here = a/(4.541kB) is the isotropic-nematic 



transition temperature in the absence of forces [19|. We see that, at temperatures below 
T n ; the molecules have spontaneous nematic order, and the force does not induce further 
order significantly. While for the molecules initially in paranematic states, the applied force 
field will induce a first-order phase transition, which means S jumps discontinuously to a 
higher value at a certain critical force fc(T*)- At even higher temperatures, nematic field is 
weaker and the effect of the force is less dramatic. Interestingly, the a — S curves at different 
temperatures are qualitatively similar to the stress-orientation curves given by a much more 



complex nematic elastomer theory [25| (Fig. 5 in Ref. [25]). 

To compare with the mechanical experiments of the SDS, we give the expressions for stress 
and strain in our theoretical framework. Apparently, the stress a of a bulk is a = Fj A = Nf: 
F is the force on the surface of the bulk, A is the area of the surface, and is the number 
of molecules per area. The strain e of the bulk is defined as e = [L(f) — L q ]/Lq, where L(f) 
is the length of the bulk along z-axis when the force field f is applied and we can take it as 
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L(f) = (l\ cos6\), and L = L(f = 0) = 1/2. Then the strain e is e = 2(| cos#|) - 1. We 
show e versus a at different temperatures in Fig. 2(c). We see that at temperatures below 
T ni , the strain grows smoothly with the stress. While for temperatures just above T ni , the 
strain grows with the stress in almost a linear way under small forces, and then a jump in 
the strain occurs at the critical force fc(T*), after which the strain increases smoothly with 
the stress again. At even higher temperatures, the jump is replaced by a smooth increase in 
strain, but there is a plateau in a certain range of force. 

In our model, the reduced temperature T* is an essential parameter, and we need to 
choose a proper value for it in order to predict the stress-strain curve of the SDS. Experiments 
showed that the solution from which the SDS was drawn was in liquid crystalline state at 



ambient temperature and pressure 



y, y, q, 



171 ]. while the orientation in the amorphous 



region of solid silk was very low 14j. Thus, we assign the isotropic-nematic transition 
temperature T ni of the cross-linked chain network slightly lower than the room temperature 
T r . Namely, the SDS is in paranematic state at ambient temperature and T* is just above 
1. The a — e curves in the paranematic states in our model indeed exhibit main features 
of the stress-strain relation of the SDS: there is a linear increase in stress with strain at 
small values, and then at a certain strain and afterwards, the material becomes softer with 
lower Young's modulus 18|, [26|. We then reveal the beginning of the isotropic-nematic phase 
transition as the yield point. We see that the curves with T* = 1.01 and 1.02 agree well 
with the measured curve in the beginning linear region and the yield point. But because 
the actual deformation process of the silk is more complicated and additional factors may 
be involved, such as the viscoelasticity, defects and poly-domain effect 20], the curve with 
T* = 1.1 agrees better with the overall stress-strain measurement topologically. We calculate 
the yield point by choosing curves with T* = 1.01 and 1.02. We get the yield strain e y ~ 0.04, 
the yield stress a y = aNk^T/l pa 8.4MPa, and the Young's modulus at the linear region 
E = a y /Ey 210MPa, given a ~ 0.2, N/l ~ 10nm" 3 , and k B T r ~ 4.1pNnm. These results 
agree with experimental data III, |26| [Fig. 2(d)] satisfactorily. 

We notice that our results agree much better with the mechanical properties of the silks 
with low spinning speed. That is because the spinning speed can induce a low orientation 
in the amorphous region which makes the silk more stiff. Since this additional order in the 
amorphous region is not taken into account in the current work, the silks in our model are 
generally a bit softer than the silks with high spinning speed. Another thing needs pointing 
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out is that we predict there is a phase transition at the yield point which is supported 
by a few experiments. For instance, in the polarized FTIR spectroscopy experiment by 
Papadopoulos et ai, the orientation of some components in the amorphous region increased 
by 0.3 when the strain reached 24%. Besides, we would like to discuss about the isotropic- 
nematic transition temperature T ni of the cross-linked chain network. If we choose T* = 1.1 
in our calculation, the transition temperature T n \ pa 270K, which is reasonable. 

In addition to describing the stress-strain relation of the SDS, our simple theory can also 
qualitatively account for the drop of the stress in the wet SDS, i.e. the supercontracted 
SDS. We take Lq and R® as the initial length and radius of the silk, and L and R as those 
under stress. Under the assumption of volume conservation we have ttR^Lq = irR 2 L, so 



R/Ro = y 1/(1 + e). The free energy of the bulk can be written as 

F = U-f ext (L-L ) + 27rRL 1 



U - f ext (L - L ) + 2nR L 0l VT+^, 



(6) 



where U is the internal energy of the bulk, f ext is the external force on the bulk and 7 is the 
surface energy coefficient. Minimizing F with respect to e, we get 

_ __f__ _ _J_?U 7 m 

When the silk is immersed in water, the surface energy coefficient 7 increases, so with 
the same stress a we get a bigger strain. Thus our theory can predict the softening of 



11 



13 



21 



2n 



29) 



supercontracted silk, an effect observed in many experiments 

In conclusion, we investigate the mechanical properties of the SDS from a point of view 
of the LC continuum theory. We found out that the deformation process is a force-induced 
isotropic-nematic phase transition process. Remarkably, such a simple model with Maier- 
Saupe theory is able to reproduce the stress-strain curve of the SDS, get the yield point, 
and qualitatively interpret the drop of the stress in the supercontracted silk. 

This work is supported by the National Innovation Research Project for Undergraduates. 
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FIG. 1: (a). A schematic diagram of the structure of the dragline silk. The bold lines represent 
the /3-sheet crystals, and the thin lines represent the polypeptide chains in the amorphous region. 
The z-axis is along the silk axis. (b). The coordinate system of the nematics. n is the director of 
the nematics, u is the director of the mesogenic molecule, and 9 is the angle between the long axis 
of the molecule and the silk axis z. 
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FIG. 2: (Color online.) (a) Comparison of a typical measured dragline silk's stress-strain curve 
(black solid line) with theoretical curves evaluated by the freely jointed chain (red dash dot line) 
and the hierarchical chain model (olive dash line) [After T. Vehoff et al.]. (b) The orientation order 
parameter S as a function of a(= fl/k^T). (c) The stress-strain curves at different temperatures 
T*(= T/T n i). (d) The stress-strain curve of the SDS spinned with the speed of lmms -1 [After N. 
Du et al.]. 
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